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Abstract 

We have completed theoretical predictions for positronium energy levels 
through order m by the calculation of the spin independent, radiative recoil 
correction. This contribution is significant and amounts to 10.64 MHz for the 
IS state. We further perform detailed comparison of theoretical predictions 
to experimental results for IS — 2S and 25" — 2P transitions. 
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Positronium is an unique hydrogenic atom to study quantum electrodynamic (QED) 
effects in bound systems. Small mass of point-like leptonic constituents: electron and 
positron, causes that the strong and weak interaction effects contribute at the negligible 
precision level. The energy spectrum may therefore be predicted with the accuracy limited 
only by the complexity in the higher order QED calculations. The equal mass of the electron 
and the positron requires a special theoretical treatments to incorporate relativistic and 
recoil effects on the same foot. Moreover, the annihilation of positronium affects also its 
energy levels. Therefore, the calculations of higher order corrections presents a challenge 
in the development of the quantum electrodynamics. Recently a significant progress has 
been achieved by the complete calculation of a single photon annihilation contribution to 
the positronium hyperfine splitting (hfs) by Adkins et al. [|l| and Hoang et al. and a pure 
recoil contribution to S-levels energies. In this letter we report on the calculation of the 
last unknown correction to positronium S-levels in the order of m. It is a spin independent 
radiative recoil contribution. Having it evaluated, we give theoretical predictions in Table |, 
for six transitions in positronium which are experimentally known with a high accuracy. In 
conclusions we emphasize that current predictions are now more accurate than experimental 
values and in some cases a disagreement is met. 

The radiative recoil corrections are in general difficult in treatment. A rigorous deriva- 
tion starts from the Bethe-Salpeter equation, and incorporates radiative corrections in the 
kernel |]15[. In the order of a{ZaYm a simplified treatment is sufficient. One finds an 
effective interaction potential between the electron and the positron from the corresponding 
S-matrix amplitude. When this amplitude is infrared divergent, than the separate treatment 
is necessary for small photon momenta. An illustrative example is the one-loop contribution 
to the hydrogen Lamb shift of order a {Z a)'^m, where the low energy part leads to the well 
known Bethe logarithm. In a higher order a {Z a)^ m this low energy part is simply absent. 
In the case of positronium we calculate the forward two-photon exchange electron-positron 
scattering amplitude at zero spatial components of the momenta. The nonperturbative ef- 



fects enter only in the lower order. This simplified treatment has been justified in Ref. [16 



Similar calculations have also been performed for the positronium hfs in Ref. |Tj], and for 
radiative recoil corrections to hydrogen Lamb shift in Ref. |T8[ . 



We start the calculation from the vacuum polarization effect. It modifies one of the 
photon propagators, what gives a small energy shift. Our result is |Ig[] : 



AE{nS) = 0.00901 , 1 

n-^ Vo6 27 TT^ / n-^ 



in agreement with the former result in Ref. [20|. The electron and the positron self-energy 



correction to positronium S-levels could be written in the form pi 



AE{nS) = --—^ / Ciq) , (2) 



nn'^ J rr'^t 



where is a self energy correction to one-particle amplitude, and it was calculated ana- 
lytically in Ref. [M. The integration is done along the Feynman contour. The expression 
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in Eq. (j^) requires subtraction of terms, that were included in the leading order self-energy 
contribution of order m, 

^^_^ -f + |ln[(g^-2mgo)/m^] _ 
TT (qr^ — 2 m go)^ (q^ + 2 m go) 

The integral in Eq. (0) with the subtracted C is finite and amounts to 

AE(nS) = — ^0.5612(1). (5) 

All other corrections up to the order m have already been calculated, so we could present 
an improved theoretical predictions. 

The main structure of the positronium spectrum is obtained from the nonrelativistic 
Hamiltonian. Leading relativistic effects of order m are given for arbitrary positronium 
states by Ferrel formula. Corrections of order appear including Lamb shift-like effects. 



recoil and annihilation contributions p2| , p3| . Final results for them can be found for example 
in Ref . . Recently Khriplovich and coworkers have calculated all m corrections to P- 
levels in Ref. |2^ using the Breit formalism. In the case of ^-states the calculation is much 
more complex. The logarithmic terms ln(a) m come from the single photon annihilation 
channel |2^, and from spin dependent part of photon exchange contributions |2^. The spin 



independent part, as it was found in Refs. [P4| , P8| , does not lead to ln(a) terms. Values 
for nonlogarithmic corrections are displayed in Table We divide them into 1-, 2-, 3- 
photon annihilation terms, and 0-, 1-, 2- radiative loop exchange terms. Additionally, 
photon exchange terms have spin independent and spin dependent part, proportional to 
operator, that leads to the hyperfine splitting. 
The photon annihilation terms have been calculated for the hfs. The contribution to 
energy levels could be found, by noting that 1- and 3-photon annihilation contribute for the 
orthopositronium and 2-photon for the parapositronium only, see Table The one-photon 
annihilation has been investigated in detail in Refs. p9|-pT| and the calculation has been 
completed recently. Most of the annihilation terms are state independent, i.e. they have 
only prefactor. The one-photon annihilation brings a nontrivial state dependence. The 
complete formula as derived by Hoang et. al in Ref. is 

m ( 1 



A£;(n3^i) = - 0.1256487 



+ ^[lnH-v[/(n) + vl/(l)] 

1-n 371-^2] 
+ ^-7^— ^ ' (6) 



24 n 96 n 

where \&(n) is the logarithmic derivation of the Euler F-function. The value for n = 1 has 
been independly calculated by Adkins et al. what forms a crucial check for this com- 
plicated calculation. Since the state dependent part is surprisingly large we recalculated 
it and got an agreement with Eq. (^. The state dependence of the vacuum polarization 
contribution was calculated also in Ref. [^. Two-photon and three-photon annihilation con- 



tributions are state independent and the final results after correcting previous calculations 
are presented in Refs. and respectively. 
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The photon exchange contributions to the hyperfine structure is presented as a spin 
dependent part in Table ||. Two-loop radiative term is given by cP' part of (1 + Og)^, where 
Oe is the electron anomalous magnetic moment. Single radiative loop exchange, i.e. radiative 
recoil correction [a iZaf' m) was found in Ref. [jl^. We have checked this result using the 
method of Ref. ||18[ and obtained agreement. Pure recoil contribution ((Za;)^m), which is 
state dependent, was evaluated in Ref. [1^ for = 1, but recently it has been recalculated 



for all nS states in Ref. [Ill] 



AE{nS) 



a m 



s_s+ \ - ln(«"^) +0.3767(17) 



+ -[ln(n)-^(n) + ^(l)] 



+ 



7 1 — n 1 — n 



12 n 



2n2 



where 



1/4 for triplet, 
-3/4 for singlet. 



(7) 



(8) 



These two results [0.167(33) versus 0.3767(17)] are in straight disagreement at n = 1. In 
Table |I| and below, we include results of Ref. [|14| , but in Table | both contradicting results 
are presented for the ground state hfs. 

The spin independent part of the photon exchange contributions can be also classified 
as a'^{Za)'^ m, a(Zaf m and iZaf' m terms. The first one, a single photon exchange, can 
be found from the well-known two-loop expression for the hydrogen Lamb shift, taking into 
account the reduced mass effects and multiplying the self-energy by factor 2. The radiative 
recoil contribution aiZaf'm is calculated in this work. It is a sum of Eqs. (^ and (0), see 
Table |I|. As well as in case of the hjs [17| the self energy contribution is relatively large 
comparing to all other terms. The last term ((Za)^m), pure recoil contribution, has been 
recently calculated in Ref. p[. This correction is state dependent and amounts 



a m 



-0.31056(63) 



n 1 
+ 



n 



69 



4n ' 3?T,2 512 
The sum of all constant terms for IS and 2S states is (see Table H) 
E^^\lS) =ma^ [0.16104(64) - 0.3868(22) s_s+] , 
Ei^\2S) = [0.34554(64) - 0.0992(22) . 



(9) 

(10) 
(11) 



The significant state dependence comes mainly from the 1-photon annihilation contribution. 
As an example the term contributes 0.75 MHz for 2^5*1 state. 

Unknown higher order terms limit the precision of theoretical predictions. The double 
logarithmic term ln^(Q;)m is known only for the hyperfine structure and annihilation 
contributions |3^ , and is included in the theoretical value for the hfs in Table |. In the case 
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of the spin independent part the double logarithmic contribution is unknown. We estimate 
it by 1 MHz for the ground state by scaling the well known non-recoil correction to the 
hydrogen Lamb shift with the reduced mass factor (/x/m)^. 

The most sensitive on the current calculations is the theoretical prediction for the IS*— 25* 
transition: 



E{2^Si) - E{l^Si) = 1233 607221.0(1.0) MHz, 



(12) 



where we use the values of the constant: a~ = 137.035 999 93(52) and cRoo = 
3 289 841960.394(27) MHz [0. The result is in a moderate agreement with the experi- 
ment P 



E{2^Si) - E{l^Si) = 1233 607216.4(3.2) MHz. 



(13) 



This new correction of the order a{Za)^ m (Eqs. (^) and (|^)) contributes significant amount 
of -9.31 MHz to this transition. The comparison of theoretical predictions to a less precise 
15 — 25* measurement as well as for the 2S — 2P experiments is presented in Table |. 
In all cases theoretical predictions are more precise. The 25* — 2P transitions are weakly 
affected by the discrepancy in theoretical predictions, i.e. by about 0.01 MHz. Experimental 
results of Refs. and |l^ are in agreement with theoretical predictions, while the results 
of Refs. 



and [ITT 



are in a slight disagreement. The hyperfine structure measurements 
are in agreement with the Caswell and Lepage result |2^ and in 3.4 standard deviation 
disagreement with Ref. [T3|. This disagreement could not be resolved by comparison of the 
IS — 2S transition, since it is not sufficiently sensitive on hfs contributions, see Eq. (|^) 

The evaluation of {Za)^m contribution was a complicated task. Both calculations of 
positronium hfs at this order: Refs. p6[ and |]14| reproduce in the limit of heavy nucleus mass 
M known results both in the orders of (Za)^m^/M and {Zafw? jM'^ Moreover, the 
result of Ref. |14] is in agreement with the state dependent terms of these orders calculated 



in Refs. |^] and |^ respectively. The discrepancy between these works will have to be 
resolved before final conclusions could be drawn. 
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TABLES 



Transition 


Exp. Refs. 


Experiment [MHz] 


Theory [MHz] 


Difference [MHz] 


23^1 - l^Si 


i 
i 


1233 607 218.9 (10.7) 
1233 607 216.4 (3.2) 


1233 607 221.0 (1.0) 


2.1(10.7)(1.0) 
4.6(3.2)(1.0) 


l^Si - l^So 


i 



203 387.5 (1.6) 
203 389.10 (0.74) 


203 388.20 (0.80) 
203 392.12 (0.50) ^ 


0.7(1.6)(0.8) " 
4.6(1.6)(0.5) 
-0.90(0.74)(0.80) " 
3.02(0.74)(0.50) 


2^Si - 2^Po 


i 


i 

101 


18504.1 (10.0)(1.7) 
18499.65 (1.20)(4.00) 


18498.42(0.13) 


-5.7(10.0)(1.7) 
-1.2(1.2)(4.0) 


2^Si - 2^Pi 


i 


13001.3 (3.9)(0.9) 
13012.42(0.67)(1.54) 


13012.58(0.13) 


11.3(3.9)(0.9) 
0.2(0.7)(1.5) 


2^ Si - 2^P2 


i 

in 


8619.6(2.7)(0.9) 
8618.4(2.8) 
8624.38(0.54)(1.40) 


8626.88(0.13) 


7.3(2.7)(0.9) 
o. 5(2.0 j 
2.5(0.5)(1.4) 


2^Si - 2^Pi 




El 

101 


11181(13) 
11180(5)(4) 


11185.54(0.13) 


5(13) 
6(5)(4) 


2^Si - 2^So 


not measured yet 


25424.69(0.06) ^ 





TABLE L Comparision of experiments to theoretical predictions. In case of the fine structure 
we give statistical and systematic error separately because some data in Table are correlated. The- 
oretical results for the ground state hfs were given in Refs. [|T| and [|2|. References to contradicting 
theoretical works are: ~ [|^]; ^ - |l^. Theoretical uncertainties are dominated by estimates 
for unknown higher order terms. The differences for 25" — 2P are presented without theoretical 
uncertainties, since they are significantly smaller than experimental ones. The levels are labeled 
as n^-^+^L/. 
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Contribution 


K 


AE{l^So) [MHz] 


AE{l^Si) [MHz] 


I7 


-0.12565 (1 + 





-2.34 


27 


0.03248 (i - s_s+) 


0.61 





37 


-0.05194 (1 + s_s+) 





-0.97 


0? (Za)'^ m 


0.02647 - 0.01374 


0.69 


0.43 


a{Za)^ m 


0.5702(1) - 0.5397(14) 


18.19(2) 


8.12(1) 


(Za)^ m 


-0.31056(63) + 0.3767(17) 


-11.07(2) 


-4.04(1) 


total 


0.16104(64) - 0.3868(22) s_s+ 


8.42(3) 


1.20(2) 



TABLE II. All constant a m contributions to IS state. The spin dependent part is similar to 
that in Ref. [||], although we use a different naming. Results are presented both, in relative units 
oi K m and in frequency units. 
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